Double charmonium production in exclusive bottomonia decays. 
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This paper is devoted to the leading twist exclusive bottomonia decays with double charmonium 
in the final state. Using models of the twist-2 charmonia distribution amplitudes the widths of these 
decays are calculated within light cone formalism. In addition, the processes under consideration 
are studied within nonrelativistic QCD. In our analysis we have found that the production of some 
of the P-wave charmonia mesons with L z 7^ is allowed already at the leading twist approximation. 
This means that the selection rules which predict the suppression of such decays are violated. The 
mechanism which lies behind this violation is discussed. 
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O I. INTRODUCTION. 

(N : 

Double charmonium production at B-factories has been one of the most challenging problem in quarkonium physics 
for many years. This problem appeared after the measurements of the cross sections of the processes e + e" 
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J/^rjc, J/^fr]' c , tp'rjct ^'v'ci J/^XcOi "4>'XcO at Belle [2, HJ an d BaBar @ collaborations, which were approximately by an 
(f-) ' order of magnitude larger than the leading order nonrelativistic QCD (NRQCD) predictions 0,[1,Q. There were many 
attempts to resolve this discrepancy 0, H, H, OH EH, [Hi [H, EH EH Es E3- Lately, it was shown that the agreement 
between theor y an d experiment can be achieved if one takes into account radiative and relativistic corrections to the 
cross sections [H, Eil H3| ■ 

Intensive study of double charmonium production at B-factories has already led to the considerable theoretical 
progress in understanding of hard exclusive processes with charmonia production. Further study of these processes 
can improve our knowledge about QCD dynamic of the exclusive processes. Moreover, since double charmonium 
production is sensitive to the parameters of charmonia wave functions, these processes can be used to study the 
structure of charmonia mesons. 

Another processes which can be used to study the dynamic of hard exclusive processes and structure of charmonia 
mesons are double charmonium production in bottomonia decays. These processes are very similar to the processes of 
double charmonia production at B-factories since the masses of bottomonia are very near to the center mass energy of 
ly-j ' e + e~ beams at B-factories M b i ~ y/s = 10 GeV. At the same time exclusive bottomonia decays have very important 
■^j- , advantage: due to the different quantum numbers of the initial bottomonia one has access to the final states whose 
f^*) ' production is suppressed at B-factories. 

. Although exclusive bottomonia decays with double charmonium in the final state are very interesting from theo- 
retical point of view, thus far only few papers were devoted to the study of some of these decays [2l|, US US HH . The 
Q\ \ present paper is devoted to the study of the leading twist double charmonia production in bottomonia decays. To 
carry out this study we are going to apply light cone formalism (LC) [IS HI]. Within this approach the amplitude of 
hard exclusive process can be separated into two parts. The first part is partons production at very small distances, 
which can be treated within perturbative QCD. The second part is the hadronization of the partons at larger distances. 
For hard exclusive processes it can be parameterized the by process independent distribution amplitudes (DA) , which 
can be considered as hadrons' wave functions at light like separation between the partons in the hadron. It should 
be noted that DAs contain information about the structure of charmonia mesons. To do the calculation of the decay 
widths we are going to apply models of the S- and P-wave charmonia DAs proposed in papers [27|, [28|, [29|, [3fJ . We 
also apply NRQCD to study the same processes. 

This paper is organized as follows. In the next section the description of the approach applied in this paper is given. 
In this section we show that the production of some P-wave charmonia mesons with L z ^ is allowed already at 
the leading twist approximation, what violates well known selection rules [2fj|. [3lj. The mechanism which lies behind 
this violation is discussed. In section III we present analytical expressions for the amplitudes of the processes under 
study obtained within LC and NRQCD. Numerical results of this paper are given in section IV. In the last section 
we summarize the results of this paper. 
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II. DESCRIPTION OF THE APPROACH. 
A. The amplitudes of the processes under study. 

In this paper double charmonium production in exclusive bottomonia decays will be considered. The presence of 
large energy scale M bb which is much greater than the masses of the final charmonia mesons M bb 3> M cg allows one to 
apply light cone formalism [2a ] (LC). Within this formalism the amplitude of the processes T is expanded in inverse 
powers of the large energy scale 1/M bb 

r= *o_ + _*i (1) 

M n M n+1 

In this paper we are going to consider the bottomonia decays which are different from zero already at the leading 
approximation in 1/M bb expansion. From here on these processes will be called the leading twist processes. Moreover, 
the calculation done in this paper will be restricted by the calculation of the first nonvanishing contribution in 1/M bb 
expansion. At this level of accuracy the amplitude of hard exclusive production of mesons Mi , M2 in bottomonia 
decay can be written in the following from 

T = J J d^ 2 H(^,^ Ml (6)0m 2 (6), (2) 

where £1, £2 ar e the fractions of the relative momenta of the whole meson carried by the quark- antiquark pair in the 
mesons Afi,Ma correspondingly, -ff(£i,£2) is the hard part of the amplitude, (f>Mi (£1): 4>m 2 (£2) are the leading twist 
distribution amplitudes (DA) of the mesons M±,M2- 
Now few comments are in order. 

1. The leading twist DAs </>Mi (£1), 0a/ 2 (£2) parameterize infinite series of the twist-2 operators. For instance, the 
leading twist DA of the pseudoscalar meson P parameterizes the operators Q757+-D+Q, n = 0, 1, 2.. as follows 1 

(P(q)\Qj 5l+ DlQ\0) = lfp{q + ) n+1 [ dx<f>(OC, (3) 

Jo 

where q is the momentum of the pseudoscalar meson P, fp is the constant which is defined as (P(q)\Qj^'y5Q\0} = 
ifpq^. One can also think of 0(£) as about the amplitude to find the quark-antiquark pair in the meson P with the 
fraction of the relative momentum of the whole meson £. It should be noted that DAs parameterize the nonperturbative 
effects in the amplitude. In Appendix A we collect the definitions of all leading twist charmonia DAs needed in this 
paper. 

From formulas @, (O one sees that amplitude §5§ resums infinite series of the twist-2 operators in the amplitude. 
If the mesons Mi , M2 are nonrelativistic mesons, formula ([2]) resums relativistic corrections to the amplitude T. 

2. The hard part of the amplitude H (£1,^2) describes small distance effects, which can be calculated within 
perturbative QCD. At the same time the DAs 0Mi (£i)> 4>m 2 (£2) parameterize nonperturbative effects, which take place 
at large distances. From this one can conclude that formula ^ separates the effects of small and large distances. 

In this paper we study the processes with bottomonia mesons in the initial states. Below the bottomonia mesons 
will be described at the leading order approximation of NRQCD. At this level of accuracy the amplitudes of P- 
wave bottomonia decays contain the contributions coming from color-octet states [32| . However, it is clear that the 
contributions of such states to the amplitude of hard exclusive decays are suppressed. This can be seen as follows. 
Besides the quark-antiquark pair, color-octet state contains one additional gluon. The amplitude to attach this gluon 
to the one of the outgoing charmonia is strongly suppressed since this gluon does not have enough energy. So, this 
gluon must be absorbed in the hard part of the amplitude. What leads to the suppression in a s and higher powers of 
relative velocity of bottomonia. For this reason, only color singlet states will be taken into the account in the present 
analysis. The calculation will be done using the technic of the projection operators [3. f33j| . 

As it was noted the hard part of the amplitude if (£1,^2) contains the small distance effects. At small distances 
the strong coupling constant a s is small, so one can expand the H (£1,^2) in a series over a s . However, presence 
of two strongly separated energy scale M bb ^> M cg , gives rise to the appearance of large logarithm log E\jMp ~ 
log(M bb / M^ c ) . This logarithm enhances the role of radiative corrections. The main contribution to amplitude @ 



1 Here the designation a+ = clq + a z for fourvector a M is used . It is also assumed that very energetic meson P moves along 2 direction. 
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comes from the leading logarithmic radiative corrections ~ (a s logE?/Mp) n . It turns out that these corrections can 
be taken into the account in formula ^ if this formula is rewritten as follows [25|, [26| 



To resum the leading logarithmic radiative corrections coming from all loops the scale /z should be taken of order of 
~ M bb . The hard part of the amplitude £2, /■*) should be calculated at the tree level approximation. At this 

level /f (£1, £2; /•*) depends on the renormalization scale \i only through the running of the strong coupling constant 
a a {p). The rest of the leading logarithms are resummed in the DAs 4>Mi (£1, h)4>m 2 (£2, M) using renormalization group 
method. It should be stressed that formula jlj exactly resums the leading logarithmic radiative corrections which 
appear in all loops. In the calculation we use /i = M bb /2. 

Commonly, to study the production of nonrelativistic mesons one uses effective theory NRQCD 32]. NRQCD deals 
with three energy scales m c 3> m c v 3> m c i> 2 , where m c is the mass of c-quark, v <C 1 is relative velocity of quark 
antiquark pair in charmonium. In the process of hard nonrelativistic meson production there appears one additional 
energy scale Eh ~ M bb which is much greater than all scales m c ,m c v,m c v 2 . Evidently, it is not possible to apply 
NRQCD at this scale. From the effective theory perspective, first, this large energy scale must be integrated out. And 
this is done through the taking into account renormalization group evolution of the DAs 4>Mi (£1 , A t )<^M 2 (£21 lA- 



In this section we are going to determine what processes of double charmonium production in bottomonia decays 
are the leading twist processes. To do this one should determine the asymptotic behaviour of the amplitudes of the 
processes under study in the limit M bb — > 00. The determination of the asymptotic behaviour was a subject of many 
papers [26l l3ll l35l . [3a . |37| . It turns out that this behavior is determined by the quantum numbers of the final hadrons 
and does not depend on the DAs of the final hadrons. Below we are going to follow paper (26|. The authors of this 
paper formulated general selection rules which can be used to specify the leading twist processes 

1. For a very energetic charmonium the c-quark helicity coincides with the projection of its spin into the direction 
of the hadron momentum. This rule is valid up to the corrections ~ k±/M bb , where k± is the transverse momentum 
of the quark in charmonium. 

2. For the charmonia states with L z 7^ (L is a quark angular momentum, the hadron is assumed to move along 
z-axis) the asymptotic behavior of the amplitude is power suppressed. 

3. The QCD interaction has a vector nature and at the quark-gluon vertex the quark helicity is conserved (up to 
the corrections ~ m c /M bb ). 

At the leading order approximation in the a s , the exclusive decays of the C-even bottomonia can be described by 
the Feynman diagrams similar to that shown in Fig. la. Typical diagram for the C-odd bottomonia exclusive decays 
is shown in Fig. lb. Applying selection rules (l)-(3) and diagrams shown in Fig. la,b one can prove that there are 
two types of the leading twist processes. The first one contains the processes in which the spin of the initial state 
J z — (z-axes is chosen in the direction of motion of final charmonia mesons). In this case the helicities of final 
mesons are Ai = A2 = 0. Any bottomonium meson can decay in this way. The second type of the processes contains 
the processes with J z = ±2 and Ai = — A2 = ±1. Evidently, only the \b2 meson can satisfy this condition. Note 
also, that one can introduce the quantity A = Ai + A2 and for all leading twist processes this quantity is zero. The 
amplitude of the process with A 7^ is suppressed as ~ (M c5 /Al bb ) A relatively to the A = process. 

For the bottomonia with J z = one can derive one more selection rule. To do this let us introduce the quantum 
number "naturalness": a = (— 1) S P, where P is the parity, S is the spin. If the naturalness is not conserved: 
^initial 7^ (&1&2) final ) there appears the antisymmetric tensor t^ vp<J in the amplitude. This tensor is contracted with 
the polarizations and momenta of the final charmonia. However, the polarization vector of energetic longitudinally 
polarized meson is proportional to its momentum, so for the mesons with the longitudinal polarizations there will be 
no enough fourvectors to get nonzero result after the contraction. From this one can conclude that for the leading 
twist processes the naturalness is conserved. 

Applying the rules that were presented above one can find all leading twist exclusive bottomonia decays. Some 
of these decays are the decays of the C-even bottomonia mesons and some are the decays of the C-odd bottomonia 
mesons. Evidently the C-odd bottomonia decays can proceed through at least tree gluons (see Fig. lb). From this 
one can conclude that the amplitudes of the C-odd bottomonia decays are suppressed as ~ oisjix ~ 0.07 in comparison 
with the C-even bottomonia decays. For this reason, we are not going to consider these decays. So, below the following 
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B. Selection rules. 
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FIG. 1: Typical diagrams for C-even (fig. a) and C-odd (fig. b) bottomonia decays into charmonia pair 



decays will be considered 

Vb, Xbi -> h c J/ip(ip'), Vciv'c) Xco(Xc2), Xd XcoiXci), 

XbO, Xb2 -> Vciv'c) Xcl, VciVc) Vc{rfc)> h c h c XcO Xc2, XcO XcO, 

Xcl Xcl, Xc2 Xc2 

Xb2 -> h c J/ip(ip'), XaXc2- (5) 



C. Violation of the selection rules. 



Let us consider a longitudinally polarized Xcl meson moving along z-axes. In this meson the spin S — 1 and orbital 
momentum of L = 1 of the quark-antiquark pair sum to J = 1 state. The Clebsch-Gordan coefficients [38[ of the 

longitudinally polarized Xcl meson ( J z = 0) are Cx"=+i 3 = _i = — s =+1 = l/^Ci'ios =0 = 0- This 

means that the longitudinal polarization can be realized only through the L z = ±1. However, the selection rules 
tell us that such states are suppressed. From this one can draw a conclusion: the processes with the production of 
the longitudinally polarized Xci meson (for instance, the XbO ~ > VcXci decay) are not the leading twist processes. 
Similar situation takes place for a transversely polarized h c meson. In this case one has S = 0, L = 1. Evidently, 
the state J z — ±1 can appear only if L z = ±1. So, from the selection rules one can conclude that the production 
of the transversely polarized h c meson is suppressed. Nevertheless, the calculation done within NRQCD shows that 
these conclusions are not correct. The processes shown in list (|5|) with the production of the transversely polarized 
h c meson and longitudinally polarized Xcl meson are the leading twist processes. 

To demonstrate that the production of the longitudinally polarized \ci meson is not suppressed, let us consider 
the process of the production of the i] c meson and {cc} pair \bo —> Vc{cc}. Assume further that {cc} pair has very 
small invariant mass M 5c <C Mi b . So, we deal with quasi-exclusive process for which one can apply LC technic. For 
instance, one can expand the amplitude of this process in and use the selection rules. Similarly to exclusive 

processes, the leading order contribution in X/M^ expansion for quasi-exclusive process appears due to the leading 
twist operator of the r\ c meson. Feynman diagrams that give contribution to this process are similar to that shown 
in Fig. la. It is not difficult to see that the leading order contribution to the amplitude in 1/Mg 6 expansion can be 
written as follows 

T = B^ x {crfuPw}, (6) 

where B^ u is proportional to the amplitude that describes the decay of the Xbo into two gluons, P — P175 is the dirac 
structure of the leading twist wave function of the r\ c meson with momentum p± . The tensor B^ v can depend only 
on the tensor structures q^q^ , q% 9i , ?2 12 , 1i Qi , {^'hb)^ 1 '■ Evidently, the leading twist contribution to the amplitude 
originates only from the last tensor structure. The others give terms proportional to masses. From this one can 
conclude that the leading twist contribution to the amplitude T is proportional to 

T - c%pi757''c ~ CP175C. (7) 

Now let us consider what cc states contribute to the amplitude T. To do this let us consider the expression for T in 
the center mass of cc pair and assume that cc is nonrelativistic pair. In this case one can expand the operator cj^j^c 
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in relative velocity of quark-antiquark pair using Foldy-Wouthuysen-Tani transformation [3£ 
the production can be written as follows 



1 



D 2 



2m? 



5»i J_^+ 



D X cr 



X + 0(v 3 ), 



The terms relevant to 



(8) 



where D is vector part of the gauge derivative, the index i is vector index, m c is the mass of c-quark, and x are 



Pauli spinor fields that create a quark and an antiquark respectively. Now it is clear that the operator tp^ 



1 



2m'i 



X 



corresponds to the L — Q, S = 0, J = state. So, one can recognize f? c meson contribution in this operator. From this 



one can conclude that the process Xbo ~* VcVc is the leading twist process. The operator ?/> + 



D x cr 



X corresponds 



to the L = 1, S = 1, J = 1 state. So, one sees that this operator creates the longitudinal Xci meson already at the 
leading twist approximation what violates rule (2) of the selection rules. Similarly, one can show that the transversely 
polarized h c meson can be produced at the leading twist approximation. 

Now, let us try to understand why the second selection rule is violated? To answer this question first let us recall 
how this selection rule originates. In LC leading twist distribution amplitude </>(a;) of nonrelativistic meson M moving 
along z-axes can be written in the following form 



#c,Q 2 ) 



k±<Q 2 



(9) 



where ipM(k±, x) is the hadron wave function of the mesons M , x is the fraction of momentum carried by quark. The 
meson M with L z ^= has the function ipM(k±,x) which contains the factor ~ exp (iL z ip). So, the integration over 
d 2 k± acts as a projector to the L z = state and L z ^ states cannot appear at the leading twist approximation. 
Similar arguments were used in papers [26l . [3l| to prove the second selection rule. 

Further, this point will be considered in more detail. To get L z ^ in addition to the exp (iL z tp) one should find 
somewhere the factor k±_ . It can be found in the matrix element of hard exclusive process or in the expansion of the 
bispinors of quark-antiquark pair in relative momentum (see formula (JHJ) ). The first source of the k± is not very 
important since it always leads to the suppression of the amplitude, so it is beyond the leading twist accuracy. 

Consider the case when the k± appears due to the expansion of quark-antiquark bispinors in relative momentum. 
Formula ([8]) proposes the example of the expansion of the operator 07^750 in relative momentum. As it was already 
noted, the first term of this expansion corresponds to the leading twist operator of the rj c meson. This can be seen 
from the fact that the leading twist operators among the operators of the fixed dimensionality are those which are 
maximally enhanced due to the Lorentz transformation from the center mass frame to the frame where the meson 
has very large energy. Evidently, the S^q in the center mass frame can be rewritten in a covariant way: S^q = p^/M , 
where p and M are the momentum and mass of the meson. From this one can conclude that due to the Lorentz 
transformation the first operator in formula ([8]) will be enhanced by the first power of large 7-factor, what exactly 
corresponds to the leading twist (see list of the leading twist DAs in Appendix A) . 

The second term in formula © corresponds to the longitudinally polarized Xci meson. This term contains gauge 
derivative which becomes tp + D±x ~ So, due to this term one has additional k± which leads to the nonzero 
contribution of the L z 7^ states. Here few comments are in order. First, note that in the bispinors the term ~ k± 
goes with Pauli matrices a j_ which flip the spin of quark or antiquark. This leads to the violation of the first selection 
rule. So, the second term leads to the violation of the first and second selection rules. Second, it should be noted 
that the states (S — 0, S z — 0) and (5 = 1, S z = ±1) cannot be produced simultaneously without one spin flip. This 
spin flip can be realized due to the violation of the first rule or the third rule. The violation of the third rule always 
put the process beyond the leading twist accuracy, since it is based on the vector nature of the quark-gluon vertex. 
So, the violation of the first selection rule in the wave function is needed in order to avoid the violation of the third 
selection rule. 

Now let us do Lorentz transformation of the second term in formula ([8]) from the center mass frame to the frame 
where the meson has very large energy. From the first sight, transverse components of fourvectors will not be 
enhanced by 7-factor. This means that the violation of the first and second selection rules leads to the appearance 
of the transverse components of fourvectors in the operators what puts the amplitude of such process beyond the 
leading twist accuracy. This conclusion is in agreement with the results of papers [26|, [3l[ . However, note that the k± 
and <7j_ appear in formula © in the form of vector product. This means that covariant expression of formula (J5J has 
the form 



c7m75C ■ 



M V 



1 



D 2 



2m? 



1 Pv ,+ 

m r M 



D a ct p 



X + 0(v 3 ). 



(10) 
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It is seen from this equation that although transverse components are not enhanced the structure e vtl \ p p v compensates 
this drawback and returns the second term to the set of the leading twist operators. 

At the end of this section we would like to note that the violation of the selection rules in the case of the transversely 
polarized h c meson is similar to that of the longitudinally polarized Xci meson. So, the conclusion of this section is 
the production of the longitudinally polarized Xci meson and transversely polarized h c meson is nonzero at the leading 
twist approximation due to the appearance of the special structures in the corresponding DAs which compensate the 
penalty for the violation of the first and second selection rules. 

D. Description of the calculation procedure. 

To calculate processes ([5]) at the leading twist approximation of LC one can apply the following rules. As it was 
noted above, the bottomonia mesons will be treated at the leading order approximation of NRQCD using the technic 
of the projection operators 0, [33|. Applying this technic one can calculate the amplitude of the decay of some 
bottomonium into two quark-antiquark pairs. The total momentum and relative momentum of the first pair are p\ 
and The total momentum and relative momentum of the second pair are pi and £,ipi- At large distances these 

pairs become M\ and M% charmonia mesons with the momenta p\ , P2 ■ 2 In LC the hadronization of quark-antiquark 
pair is described by DAs. To calculate the amplitude of meson production one should replace bispinors of quark- 
antiquark pair vu in the amplitude by the corresponding distribution amplitude. The list of the leading twist DAs 
needed in the calculation can be found in Appendix A. After taking the traces over dirac and color indexes one can 
get analytical expression for the amplitude. 

Another independent approach that can be used to calculate the amplitude of double charmonium production in 
bottomonia decays is NRQCD. In this approach the final charmonia mesons are treated as nonrelativistic states. In 
this paper we are going to apply NRQCD at the leading order approximation in relative velocity in charmonia. It 
should be noted that this approximation cannot be considered as reliable. This conclusion can be drawn from the 
experience obtained in the study of double charmonium production at B-factories |20| , which tells us that relativistic 
and radiative corrections to double charmonia production can be very large. Nevertheless, in this paper the leading 
order approximation of NRQCD will be used to get independent estimation of the widths of the exclusive bottomonia 
decays under study. 

In addition to the estimation of the widths of the bottomonia decays, NRQCD can be used to study very important 
question of LC: the estimation of the corrections to the leading twist approximation (power corrections). To estimate 
power corrections one can apply the idea of duality of NRQCD and LC descriptions of the hard exclusive nonrelativistic 
mesons production: if the amplitudes of the hard process under study obtained within NRQCD and LC are expanded 
in powers of inverse hard energy scales and relative velocity one will get series equal to each other. It should be noted 
that these is no strict proof of this statement. Moreover, the general proof for the NRQCD factorization of exclusive 
process involving P-wave heavy quarkonium is not available (see, e.g., [34j]L There may be large soft contributions at 
order m c /mb to these decays. However, one can expect that this statement is true since the amplitude in NRQCD 
and LC can be expended in series of equivalent operators. Assuming that both theories can describe experiment one 
can expect that these expansions in both theories coincide. Note also that we have checked that all expressions for 
the amplitudes under study obtained in LC and expanded in relative velocity coincide with that obtained in NRQCD 
and expanded in l/M^g. 

Now it is clear how it is possible to estimate the size of power corrections in LC. To do this one should take the 
leading order NRQCD prediction for the amplitude and expand it in M cS /M b i ~ m c /mf,. The first term of this 
expansion can be reproduced by the leading twist approximation of LC expanded in relative velocity in charmonia. 
The second term in this expansion can be reproduced by the power corrections to the leading twist contribution. 
So, the second term of NRQCD expansion in ~ m c /mb gives the estimation of power correction to the leading twist 
approximation. Below this approach will be used to estimate the error of the calculation. 

III. ANALYTICAL EXPRESSIONS FOR THE MATRIX ELEMENTS. 

In this section analytical expressions for the amplitudes and the widths of the double charmonium production 
in bottomonia decays obtained within the leading order approximation of LC and NRQCD are given. Before the 
analytical expressions are given, let us introduce some useful designations. First, we introduce the definitions of the 



2 It should be noted that at the leading twist approximation one can disregard the masses of quarks and hadrons. 
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constants of the vector (V — J/^,ip') and pseudoscalar (P = r] c ,ri' c ) mesons that parameterize the following matrix 
elements 



(F(p,e)|O 7a C|0) = ttM v e a 
(V(p,e)\Ca a pC\0) p = fy(ii)(p a ep - p^e a ) 
(P( P )\C lalb C\Q) = if£p a . 

Next, we introduce the following constants for the P-wave charmonia mesons Xco, Xdi Xc2i h c 



(11) 



(Xco(p)\C la (-iD„)C\0) 
(Xci(p,e)\C lal5 C\0) 
(h c (p, e)\C lal5 (-iD„)C\Q) 
(Xc2(p,e)\C la (-iD„)C\0) = fy(ji)M. 



fxi M x.i e a, (xci(p,e)\C<x a p(-iD v )C\0) = /Ji (fi)e a p p xe p p x p^, 
fk(n)M hcPa e„, (h c (p,e)\Ca a0 C\O) = .g (^)e a0pX e p p\ 



(Xc2(p, e)\Ca a /3(-iD v )C\0) = if^ 2 (fi)(p a ep v - P0€ au ) (12) 



It should be noted here that for the mesons Xci,h c the polarization e is described by the four vector e M , for the 
Xc2 meson the polarization e is described by the tensor e^ v . The superscript T in formulas (|12[) means that 

the corresponding meson is transversely polarized (A = ±1). The superscript L implies that the meson has helicity 
A = 0. Except the constants fy,fp,fi\, all constants in formulas ([HI) . O are scale (fi) dependent quantities. The 
anomalous dimensions of these constants can be found in papers 

In addition to constants (|TT|) . (fT2l . one needs the bottomonia NRQCD matrix elements (0| fc ), (Op) which are 
defined in [32]. The matrix element (OSj*) is proportional to the S-wave bottomonium radial wave function at the 
origin. The matrix element (Op ) is proportional to the of the derivative of P-wave bottomonium radial wave function 
at the origin. 

The decays of the <rj b meson. The decays r\ h -> h c J/ip(i//), rj c (i c ) Xco(Xc2), Xd Xco(Xc2) are the leading 
twist decays of the rjb meson. According to the discussion in the previous section at the leading twist approximation 
of LC only the production of the longitudinally polarized mesons is allowed. The same is true for the decays of all 
bottomonia except the Xb2 meson. Corresponding amplitude for all the leading twist decays of the rjb meson can be 
written as 



M(rj b -> MiM 2 ) 



27 V m b m l ' 



(13) 



where /f" 2 are the mesonic constants of longitudinally polarized charmonia and 

l 



6+6 



(i-6 2 )(i-6 2 )(i + 66) 



0f(6)02(6), 



where <f)f (£i)(f>2 (6) are the final charmonia DAs. It is not difficult to get the amplitude of any leading twist decay 
of the rjt, meson using formula (| 13|) . For instance, to get the amplitude of the process rjb — > VcXc2 one should take 
fi = fp,f% = /v2 an( l longitudinal DAs (see Appendix A) of the corresponding mesons. 

Using this matrix element it is easy to obtain the expression for the rjb — > M\M 2 decay width within light cone 
formalism: 



r (r]b — > M1M2) = 



5127r 3 a4 (Of 



^3„,4 lnbb\ r fL fh 

J I J 2 Arib) 



729 



n 

,2 2 o 



Within NRQCD one can get the following expressions for the width of the processes under study 

, I 2 

F( Vb -> M X M 2 ), 



5l2Tr 3 a 4 s (Of* 



T(ijb — > M\M 2 ) = 



729 



hit, 



fl 



NRQCD f NRQCD 
J2 



hit, 



where 



/•NRQCD _ f NRQCD _ 

Ji — Js — 



'{0% 



(14) 



(15) 
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for S*-wave charmonium states and 



fi 



NRQCD 



j-NRQCD 
JP 



'(Of) 



mi 



(16) 



for P-wave charmonium states. The factor F[r] b — > MiM 2 ] for different final states is 

F[i]b — > h c ip] = l + 32r 2 , 

. 1 16r 64r 2 

2 16r 32r 2 



f[f?6 -» XcOXcl] 

F[rj h Xc2Xd} 



3 3 3 ' 

2 32r 40r 2 800r 3 

3 + "3" + ~~3 ~ ' 

4 32r 104r 2 160r 3 
3 ~ IT + _ 3 ~ 



where r — m c /mt,. 

The decays of the \bo meson. The decays xw -» Vc(v' c )Xci, Vc(v' c ) Vc(v' c ), J/i>(i>') J/^W), h c h Cl XcO Xc2, XmXm, 
Xd Xcii Xc2 Xc2 are the leading twist decays of the Xbo meson. The amplitude for these decays can be written as 



M{ X 



60 



MiM 2 ) 



5127r 2 a 2 
27^ \ 



°i 6 ) f L f L r . 

IP Jl .12 T (Xbo) 



rrv 



2 J o ' 



where 



r(Xbo) 



4 + ^+666+^1 

4(i-^)(i-el)(i + €i6) 2 



>i(£i)02(6O- 



The width of the Xbo - * MiM 2 decay is 



60 



MiM 2 ) = 



8192 A 4 (Q b p b ) 
2187 m 4 



A /2 AXbo) 
™2 J 



(17) 



One should remember, that in the case of the identical final mesons the width should be divided by 2. Within NRQCD 
one can get the following expressions for the width 



where 



r(x 



60 



M X M 2 ) = 



4096tt 3 c4 (of) 



2187 



t NRQCD f NRQCD 
J I J 2 



F{xbo -» MiM 2 ) 



^[X60 -» ?7cXcl] 
^[X60 -» XcOXc2] 

^[X60 -» ??c»7c] 
^[X&0 -» V^] 

^[XbO -» ftc^c] 
^[XbO -» XcOXco] 
^[X60 -> XclXcl] 
^[X60 -» Xc2Xc2] 



4 - 16r, 

1 28r 92r 2 1120r 3 
9 ~ "9" + ~~3 9 

1 + 4r + 4r 2 , 
l-4r + 12r 2 , 

i - 10r 2 - 32r 3 + 272r 4 , 

128r 3 



1600r 4 



16r 



36 



28r^ 



16r 4 



9 



537r 2 

4 - 56r + — 260r 3 + 72r 4 , 

1 28r 179r 2 340r 3 1480r 4 
9 ~ "9" + "6 9" + 9 
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The decays of the Xbi meson. The decays Xbi -» h c J/tp(ip'), T) c (rf c ) Xco(Xc2), Xa Xca{Xc2) are the leading 
twist decays of the Xbi meson. The amplitude for the processes equals 



where 



128\/27r 2 Q;2 
M(xbi ^M,M 2 ) = v — ^ 



Of ) r L fL 

IP J I J 2 Axbl ) 
S 2 0' 

ml m 2 b 



4 Xbl) = I #id& 



<^ (6)02 (6)- 



The width of the Xbi — » M\Mi decay is 



r(xfti -» AfiM 2 ) 



10247r 3 at 



2187 



(Xbi) 




Within NRQCD one can get the following expressions for the width 



IXxbi -» MiM 2 ) 



where 



2187 

-» h c tp] = 

F[xbi -> VcXco] = 

F[Xbi -» ?7cXc2] = 

^[Xfci XcoXd] = 

^[X&l Xc2Xcl] = 



./l 



NRQCD f NRQCD 
32 



mt 



F( X bi -» MiM 2 ) 



l + 4r-32r 2 , 

1 4r 
3 ~ ~3~' 

2 , 4r Iff 2 

124r . „. 2 



3 +24 ^' 
^ - 8r 2 + 192r 3 



The decays of the Xfc2 meson. The decays x&2 -> »? c (»7c) Xd, »? c (»7c) 'fcfac), J/ip(ip') J/ip(ip'), h c h c , XcoXc2, 
Xco Xco, XciXci, Xc2 Xc2, ft c J/tPW), Xa Xc2 arc the leading twist decays of the X62 meson. The decay of this 
meson is more complicated than the decays described above. The point is that large spin of the Xc2 meson opens 
the possibility to produce transversely polarized charmonia (Ai = — A2 = ±1) at the leading twist approximation. 
So, contrary to the bottomonia decays discussed above in some decays of the Xc2 meson there can be two different 
polarizations of the final charmonia: Ai = A 2 = or Ai = — A 2 = ±1. It should be noted here that it depends on the 
process what possibilities are realized at the leading twist. For instance, for the decay Xc2 — > VcXa only Ai = A 2 = 
is allowed, for the decay Xc2 — » J/^J/^f both possibilities Ai = A2 = 0, Ai = — A2 = ±1 are allowed. The state 
Ai = A 2 = is forbidden for the decay Xb2 — ► h c J/ip since in this case naturalness is not conserved. So, the only 
allowed possibility is Ai = — A2 = ±1. 

The amplitude for the decays of the x&2 meson into pair of the longitudinally polarized mesons (Ai = A 2 = 0) equals 



„w„wn 256V27r 2 a 2 



O?) f L f L . . 
I P Jl J2 AXbl) 



where 



r(X6 2 ) 
J 



1 

/ 



2(i-??)(i-e 2 2 )(i + 66) 



,2 J o > 



2^ (6)0 2 L (6). 



If both final mesons can have nonzero helicity Ai = — A2 = ±1, one should take into account the decays into 
transversely polarized particles. The amplitude of these decays is 



M{xb2^M?Mj) =- 



512ttV 



27 \ 



Of) f T f T . . 

I P Jl J2 T {Xb2) 



ml ml 2 
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where 



Axb2) 



The width of the Xb2 — * M X M 2 decay is 



= / #1*2 



^(6)^(6)- 



T( X b2^M 1 M 2 ) = 



40967r 3 a 4 (Off) 
10935 m 4 



flf2 j{Xb2) 
2 



12 



X fT fT 
Jl J2 AXb2) 
2~~ 1 2 



rnr 



Within NRQCD one can get the following expression for the width of the process \b2 —> M\M 2 

2048vT 3 a 4 (OS?) 



where 



10935 m 4 



/•NRQCD f NRQCD 
Jl 12 



F(Xb2 -> MiM 2 ) 



-F[X&2 -> ?7cXd] 

^[Xb2 -» XcOXc2] 

^[Xh2 -> ^c] 
^[Xb2 -> V^] 
-F[Xfc2 -> /l c /lc] 

^[Xfc2 -> XcOXco] 

-F[Xb2 -> XdXcl] 

^[Xfc2 -> Xc2Xc2] 

F[Xb2 -> /lc^] 
^[X62 -> XclXc2] 



4 + 8r-96r 2 , 
16 872r 208r 2 544r 3 
~9 + ~9~ + ~~ 3 IT" 



4288r 4 
IT" 



9 

= 1 - 8r + 16r 2 , 
= 13 + 56r + 48r 2 , 
= 16 - 132r + 488r 2 - 944r 3 + 800r 4 , 
4 16r 16r 2 128r 3 256r 4 
9 



9 3 9 9 

7 - 44r - 30r 2 + 340r 3 + 264r 4 , 



43 44r , 6212r 3 

= 1 286r 2 H 

9 9 9 

= 24- 72r-96r 2 , 

= 6 + 4r 2 - 488r 3 + 1072r 4 



5032r 4 
IT" 



(18) 



IV. NUMERICAL RESULTS. 



In order to obtain numerical results from the presented above analytical expressions the following numerical pa- 
rameters are needed. 

In this paper we are going to use the models of the charmonia DAs proposed in papers [U [H, H^, [3(| ■ For the 
strong coupling constant we use one-loop expression 



a a (M) 



47T 



6 1u(m 2 /A 2 



QGD) 



where b = 25/3 and A QCD = 0.2 GeV. 

In the calculation the following values of the constants f i ' T defined in equations Ijlip. (fTSj) will be used 



f L 


= 0.373 ±0.064 GeV, 






fj/i> 


= 0.416 ± 0.005 GeV, 




= 0.379 ±0.021 GeV 


/jt(2S) 


= 0.261 ± 0.077 GeV, 






/■0(2S) 


= 0.303 ± 0.003 GeV, 


/J(2S)( M .//* 


) = 0.261 ±0.042 GeV, 




= 0.093 ±0.017 GeV, 








= 0.160 ±0.015 GeV, 


fl(Mj/*) = 


0.179 ± 0.032 GeV 


f L 


= 0.272 ±0.048 GeV, 




= 0.111 ±0.020 GeV, 


fxJMj/*) 


= 0.131 ± 0.023 GeV, 




= 0.157 ±0.028 GeV. 



(19) 
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The values of the constants /f/,,,, /^ra'Ji can ^ e extracted from the leptonic decay widths of the J ftp and ip(2S) mesons. 



' J/i/>' •'V(2S) 

The values of the constants f^ c , f^ 2 S) 



were calculated in paper [2C 



The values of the constants f T 



' •' 0(2f 



can 



be found in paper [40j • The values of the constants of the P-wave charmonia mesons can be found in paper [3(| . It 
should be noted that the constants fjt^ > f\co >fh' > ^xli > -^Xc'J depend on the renormalization scale. As it is seen 
from formulas (fT9|) these constants are defined at the scale /i = Mj/^,. 
which govern the evolution, can be found in papers [30L l40l ] . 



The anomalous dimensions of these constants, 



The values of NRQCD matrix elements 
function R(r) as follows 



bb 



O b i^ can be expressed through the bottomonia radial wave 



of), 



^|fis(0)l 2 , 
Ztt 



01 



bb\ 



ll 



f \Rp(0)\ 2 - 



(20) 



In this paper the values of the |Pg(0)| 2 , |i?p(0)| 2 will be determined from the Buchmuller-Tye j43,|43| potential model. 
Thus one gets 



of) ! 



3.1GeV J 



of) i 



2.0GeV J 



In the forthcoming analysis we are not going to take into the account the uncertainties in the values of the NRQCD 
matrix elements for bottomonia mesons since these uncertainties are not very important. 

In the calculation we also need the values of constants (|15[) . ([To]) for the charmonia mesons. At the leading order 
approximation in a s and relative velocity these constants can be determined from the leptonic decay widths 



e ' e = 



4<? 2 7ra 2 r -NRQCD"! 2 



3M V 



if: 



Using experimental results T(J/^f 



5.55 KeV,r(^(2S') -> e+e") = 2.48 KeV [H one gets 



[/N RQ ° D ] 2 = 0.17 ± 0.06 GeV 2 , [/ 2 N / QCD ] 2 - 0.09 ± 0.05 GeV 2 



(21) 



(22) 



To determine the constant y^J R Q CD one can use ^ ne decay width \c0 



77 



r(xc 



77) = ^^[/r QCD ] 2 . 



(23) 



Using experimental results Br(\cO — * 77) = 2.35 x 10 4 [38| one gets 



(24) 



were calculated as follows. 



In NRQCD there are relativistic and radiative 

2\ 



The uncertainties in values 

corrections to formulas (|2"Tj) and ([23]) . The relativistic corrections can be estimated as (i> 2 )is = 0.21 [27], (f 2 )2S = 0.54 
[29l | . (v 2 )p = 0.3 [33|. The radiative corrections can be estimated as ~ a s (Mjm) = 0.25. Adding these uncertainties 
in quadrature one estimates the errors of the calculation. 

The last parameters needed for calculation of the bottomonia decay widths are the pole masses of the 6-quark and 
c-quarks. For the c-quark we take m c = 1.4±0.2 GeV. For the mass of 6-quark we take the value m& = 4.8±0.1 GeV. 

In Table IIVI we present numerical results for the widths of exclusive bottomonia decays into pair of charmonia 
mesons. In the second and third columns of this table the results of NRQCD and light cone formalism arc shown. In 
the fourth column we present the branching" fractions of the considered decays. To estimate these fractions we use 
the following expressions for the total widths of bottomonia mesons [43[ : 



bb\ 



61 



4^ 2 (Of 



9 

3C F 

7T 

N c 
AC F 

In/ 

N c 
1.6 MeV 



11 MeV, 



, (Q b P b ) 

*s 4 

2 <^> 



n f 2 
— it at 

3 
n 



_l, 
3 ' 

317 



(Os) 
s ml 



54 288 1 



0.8 MeV, 



0.2 MeV, 



16 4, A \ 

In n f 

27 9 2m b l ! 



(O 



n f 2 



(o 8 



(25) 
(26) 
(27) 

(28) 
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0.5 1.0 1.5 2.0 

FIG. 2: The width of the Xbo ~* XcOXc2 decay as a function of the mass of c-quark m c . 

where n/ = 4 is number of active flavors, A = 0.2 MeV, and (Og) w 0.0021 (Op) is the color octet matrix element 
for the P-wave bottomonia mesons. Since experimentally charmonia mesons are observed in their decays into J/ijj 
meson, it is interesting to know the widths of such processes. It is clear that they are equal to 

T[(bb) -> (cc)i(cc) 2 -> J/^J/i/j + X] = T[(bb) -> (cc)i(cc) 2 ]Br[(cc)i -> J/V> + X]Br[(cc) 2 -> J/V' + X] 

These values are shown in the last column of Table IIVI 

Now let us discuss the uncertainties of the calculation. Before we discuss how the uncertainties of the NRQCD 
prediction can be estimated, one should recall the experience gained from double charmonium production at B- 
factories. In this case the leadingorder NRQCD predictions [1,H, @ are approximately by an order of magnitude 
less than experimental results Note also that measured values of the cross sections are much larger than the 

leading order NRQCD predictions even if one takes into account the possible uncertainties of the approach From 
this fact one can conclude that it is rather difficult to calculate the uncertainties of NRQCD. So, the uncertainties 
calculated in this paper can be considered only as a very rough estimation of the real uncertainties. 

Note also that in some decays ( see section III), the amplitudes contain polynomials in r with alternating signs of 
the coefficients of these polynomials. One can expect that the uncertainty of the calculation of such decays can be 
very large. The decay Xbo —> XcoXc2 can be considered as a dramatic demonstration of this point. The width of this 
decay as a function of the mass of c-quark is shown in Fig. 2. It is seen from this figure that the width \bo - > XcoXc2 
has zero and minimum at m c ~ 1.56 GeV which is very near to the pole mass of c-quark. Evidently, the uncertainty 
in the width of this decay can be very large. 

As the uncertainty of the leading order NRQCD prediction we take the uncertainty which originates from the 
following sources: uncertainty in the pole masses of c and b quarks (the first error in the second column of Tab. HVp . 
uncertainty due to the values of constants (f2"2")l . (|24p (the second error in Tab. IIV|) . uncertainty due to the unknown 
radiative corrections, which can be estimated as ~ a s (M bi ) log m 2 b /m 2 c ~ 50% (the third error in Tab. HVp . The last 
uncertainty (the forth error in Tab. IIVI) originates from the uncertainty in the value of the A = 200 ± 40 MeV what 
corresponds to a s (m T ) = 0.34 ± 0.03 [38l] . 

The uncertainties of the results obtained within LC can be divided into the following groups: 

1. The uncertainty in the models of the distribution amplitudes, which can be estimated through the variation of 
the parameters of these models (see papers [U [H, H^, |30|) (the first error in the third column of Tab. HV|) . 

2. The uncertainty in the values of constants (fT51) (the second error in the third column of Tab. HVp 

3. The uncertainty due to the power corrections. This source of uncertainty is very important and for many processes 
this is the main source of the uncertainty. To estimate this source of the uncertainty we expand of the leading order 
NRQCD results in the ratio r = ml jm\. The first term in this expansion is the leading twist contribution which is 
reproduced within LC. We take the next nonvanishing term in the r expansion as the estimation of the size of power 
corrections, (the third error in Tab. lIVp 

Applying this approach for the estimation of power corrections and looking to the leading order NRQCD results 
(section III) one can separate all processes into three groups. The first group contains the processes (for instance, 
the decays r/b — > h c xJj,Xbo ~ > VcXdt ) for which power corrections, most probably, will not change the LC predictions 
dramatically. One can expect that for this group of the decays LC predictions are reliable. The second group contain 
the processes for which power corrections are of order of ~ 100% (for instance, r\b — > rj c x c o, Xbo —> XciXci )■ For such 
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MiM 2 


TNRQCD, eV 


Tlc, eV 


BlLC, % 


BrLcO/0, % 


ri b — > h c J/ip 


73.t^; ±38. ±37. ±29. 


520. ± 63. ± 140. ± 210. 


0.0047 


— 




Vb h c ip(2S) 


39 +6 8 ± 26. ± 19. ± 15. 


300. ± 36. ± 78. ± 120. 


0.0027 


— 




rib -> VcXcO 


56.±\l\ ± 29. ± 28. ± 22. 


49. ± 25. ± 71. ± 20. 


4.5 x 10~ 4 


— 




rib -> r/ c (2S)xcO 


30.±? | ± 20. ± 15. ± 12. 


26. ± 18. ± 37. ± 10. 


2.3 x 10~ 4 


— 




Vb -> VcXc2 


17.±|;1 ± 8.9 ± 8.6 ± 6.9 


97. ±48. ±71. ±39. 


8.9 x 10~ 4 


— 




r/b -> r/ c (2S)xc2 


9.lla;?±6.1±4.6±3.6 


51. ±35. ±37. ±20. 


4.6 x 10~ 4 


— 




Vb -» XcoXci 


11+2; ±8.4 ±5.5 ±4.4 


26. ± 13. ± 38. ± 10. 


2.4 x 10" 4 


1.1 x 10" 


-6 


% XclXc2 


4.7±i i ±3.6 ±2.4 ±1.9 


51. ±26. ±37. ±21. 


4.7 x 10" 4 


3.4 x 10" 


-5 


Xbo r?cXd 


12 .+j; 4 ±6.1 ±5.8 ±4.7 


53. ±26. ±17. ±21. 


0.0064 


— 




Xbo ->■ ric(2S)xd 


6.2±^;| ± 4.2 ± 3.1 ± 2.5 


39. ± 27. ± 13. ± 16. 


0.0048 






XbO —> XcOXc2 


0.00059±2:ooo4i ± 0.00045 ± 0.00029 ± 0.00024 


0.81 ±0.42 ±1.8 ±0.33 


10. x 10~ 5 


2.6 x 10" 


-7 


Xbo -> r? c J7 c 


49 +"3 ± 35. ± 25. ± 20. 


51. ± 25. ± 16. ± 20. 


0.0062 


— 




Xbo -»■ ricilci2S) 


26+H ± 17 - ± 13 - ± 10. 


76. ± 52. ± 24. ± 30. 


0.0093 


— 




Xbo ?7c(25)?7 c (2S) 


14 .+J; 5 ± 15. ± 6.9 ± 5.5 


30. ± 25. ± 9.6 ± 12. 


0.0036 


— 




Xbo -> J/tpJ/ip 


27.±1; 2 ± 19. ± 13. ± 11. 


79. ±3.1 ±25. ±31. 


0.0096 


0.0096 




Xbo -> J/W(2S) 


14 +|| ±9.3 ±7.1 ±5.7 


130. ±9. ±41. ±51. 


0.016 


0.0089 




Xbo - V(2S)V(2S) 


7.5+^; 4 ± 8.3 ± 3.8 ± 3. 


54. ±6.7 ± 17. ±22. 


0.0066 


0.0022 




XbO — > ftcfec 


0.094±g;g|| ± 0.072 ± 0.047 ± 0.038 


15. ±2.9 ±3.9 ±6. 


0.0018 


— 




XbO — * XcOXcO 


0.048lg^34 ± 0.037 ± 0.024 ± 0.019 


0.21 ±0.11 ±1.1 ±0.082 


2.5 x 10~ 5 


4.2 x 10" 


-9 


XbO -» XclXcl 


0.56±g;f| ± 0.43 ± 0.28 ± 0.22 


14. ±6.8 ±15. ±5.5 


0.0017 


2.1 x 10" 


-4 


t Xi>0 — ► Xc2Xc2 


0.026to;°?3 ± 0.02 ± 0.013 ± 0.01 


0.81 ±0.41 ± 1.8 ±0.32 


9.9 x 10" 5 


4. x 10" 


6 


Xbi h c J/ip 


1.2+%™ ± 0.64 ± 0.61 ± 0.49 


14. ±1.8 ±4.6 ±5.7 


8.9 x 10" 4 


— 




Xbi -» /i c -0(25) 


0.65±g;i 2 ± 0.44 ± 0.32 ± 0.26 


15. ±2.6 ±4.7 ±5.9 


9.2 x 10~ 4 


— 




Xbl -» ??cXcO 


0.24tg;^ 1 ± 0.13 ± 0.12 ± 0.097 


1.3 ±0.68 ±0.43 ±0.53 


8.3 x 10~ 5 


— 




Xbi -> ?? c (25)xc0 


0-13t ( o'.027 ± 0.087 ± 0.064 ± 0.052 


1.3 ±0.91 ±0.41 ±0.51 


8. x 10~ 5 


— 




Xbl -» ??cXc2 


°- 74 io'.i2 ± o- 38 ± °- 37 ± °- 29 


2.7 ± 1.3 ±0.42 ± 1.1 


1.7 x 10~ 4 


— 




Xbi r; c (25)x C 2 


0.39to.o63 ± 0- 26 ± 0- 19 ± 0- 16 


2.5 ± 1.8 ± 0.41 ±1. 


1.6 x 10~ 4 


— 




Xbl — > XcOXcl 


0-6lo.?s ± 0- 45 ± 0- 3 ± O- 24 


0.68 ±0.35 ±3.4 ±0.27 


4.3 x 10~ 5 


2. x 10" 


7 


Xbl -» XclXc2 


0- 13 -o.o?9 ± 0- 098 ± 0- 064 ± 0- 051 


1.4 ±0.68 ±0.43 ±0.54 


8.4 x 10" 5 


6.1 x 10" 


-6 


Xb2 -> r?cXd 


i-S-o'Jg ± O- 92 ± 0.88 ± 0.71 


17. ±8.2 ±2.7 ±6.7 


0.0069 


— 




Xb2 r? c (2S')x c i 


0.931o:1b ± 0.63 ± 0.47 ± 0.37 


11. ±7.7 ±1.8 ±4.4 


0.0046 






Xb2 -» XcOXc2 


0.58lf,;?7 ± 0.44 ± 0.29 ± 0.23 


1.4 ±0.72 ±6.1 ±0.57 


5.8 x 10~ 4 


1.5 x 10" 


-6 


Xb2 -> r? c J7 c 


1.6ioi 9 ± 1.1 ±0.78 ±0.62 


16. ±7.8 ±10. ±6.4 


0.0066 


— 




Xb2 -> ri c (2S)ri c 


0.83±g;|| ± 0.54 ± 0.41 ± 0.33 


21. ±15. ±14. ±8.5 


0.0088 


— 




Xb2 - 7? C (2S)J7 C (25) 


0.44lo 14 ± 0.49 ± 0.22 ± 0.17 


8. ±6.7 ±5.1 ±3.2 


0.0033 


— 




Xb2 -» J/lpJ/tp 


65.±1|; ±46. ±32. ±26. 


270. ±41. ±93. ±110. 


0.11 


0.11 




Xb2 - </>(2S)J/</> 


34+™ ±23. ± 17. ± 14. 


380. ± 120. ± 130. ± 150. 


0.16 


0.09 




Xb2 - V(2S)V(2S) 


18 +3.3 ±20. ±9.1 ±7.3 


140. ± 63. ± 50. ± 58. 


0.059 


0.02 




Xb2 — > ftc^c 


0.42±g;i 1 ± 0.32 ± 0.21 ± 0.17 


37. ±9.5 ±24. ± 15. 


0.015 


— 




Xb2 — ► XcOXcO 


0.014±g;gg|| ± 0.011 ± 0.0072 ± 0.0058 


0.36 ±0.19 ±0.11 ±0.14 


1.5 x 10~ 4 


2.5 x 10" 


-8 


Xb2 -> XclXcl 


0.18t°,:oB9 ± 0.14 ± 0.089 ± 0.071 


8. ±4.8 ±4. ±3.2 


0.0033 


4.2 x 10" 


-4 


Xb2 — ► Xc2Xc2 


0-2±g;gl| ± 0.15 ± 0.098 ± 0.078 


9.8 ±7.8 ±0.8 ±3.9 


0.004 


1.6 x 10" 


-4 


Xb2 -> ftcJ/V' 


7.6±i'J ± 4. ± 3.8 ± 3. 


100. ± 42. ± 24. ± 41. 


0.042 






Xb2 -» /i c -0(2S) 


4+°;™ ±2.7 ±2. ±1.6 


71. ±36. ±17. ±28. 


0.029 






Xb2 -» XclXc2 


0.32+^g ± 0.24 ± 0.16 ± 0.13 


11. ±9.6 ±0.047 ±4.4 


0.0045 


3.3 x 10" 


-4 







TABLE I: The widths and branching fractions of the exclusive bottomonia decays into pair of charmonium mesons. In the 
second column the NRQCD predictions are presented. In the third and fourth the widths and branching fractions of the 
exclusive bottomonia decays in LC formalism are shown. The last column contains the branching fractions of inclusive J/ip- 
pair production through intermediate charmonium states. The symbol "X" in this column means, that this decay is forbidden 
(for example, decay rib — ► VcXcO — > J/ipJ/ip + X is absent since 7? c meson cannot decay into J/tp). The symbol " — " (for 
example, in the case r/ b — > h c J/tp) means that the branching fraction of the corresponding decay (h c — ► J/tp ± X in this case) 
is unknown. 
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processes our results are valid up to the factor of ~ 2. The last group of processes are the processes for which power 
corrections are large (for instance, Xbo ~ * XcoXcOi Xbi ~ * XcoXci)- For these processes we can guess only the order of 
magnitude of the widths within LC. 

4. The uncertainty due to the radiative corrections. The main part of the radiative corrections to the amplitude - 
the leading logarithmic radiative corrections have been resummed within LC. This fact allows us to estimate the rest 
of the radiative corrections as ~ a s (E) ~ 20%. This is very small uncertainty, so we don't show it in Tab. IIVI 

5. The uncertainty due to the variation of A = 200 ± 40 MeV. (the forth error in Tab. HV[) . It should be noted 
also that in the calculation we took the scale of factorization /i = M bb /2. However, one can take any scale /i ~ M bb . 
The dependence of the branching ratios on the /x are rather strong since Br ~ a 2 s . For instance, if we change fx from 
M bb /2 to M bb the branching ratios will be changed by 40 — 50 %. 

V. CONCLUSION. 

In this paper the leading twist double charmonium production in exclusive bottomonia decays was considered. The 
decays of the C-odd bottomonia are suppressed by the factor ~ a s /ir in comparison to the decays of the C-even 
bottomonia. For this reason we considered only the leading twist decays of the C-even bottomonia. Applying light 
cone formalism with the models of the leading twist charmonia distribution amplitudes (27l . I2H [2^ . |30| we calculated 
the amplitudes and the widths of the corresponding processes. In addition, we calculated the widths within the 
leading order NRQCD. 

During the calculation we found that the production of the longitudinally polarized Xci meson and transversely 
polarized h c meson with L z 7^ is nonzero already at the leading twist approximation. This fact tells us that the 
second selection rule ( see section II ), which predicts the suppression of such processes, is violated. We considered 
the mechanism which lies behind this violation and found that this violation results from the rather special Lorentz 
structure of the corresponding distribution amplitudes. 
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APPENDIX A: DISTRIBUTION AMPLITUDES. 



The leading twist distribution amplitudes needed in the calculation can be defined as follows: 
for the pseudoscalar mesons P = rj c ,rj' c : 



fp ^ 



P(p) Ql(z)[ Z ,-z}Q^(-z) 0) = (p75)/3a 4 
for the vector mesons V — J / ^ , i/j' '• 



^ (pz) 0p(C;m), 



(v(p,e, =0 )\Q* a (z){z,-z}Qi p (~z)\o) = {p) ti J-ff j dee*<*»<ft{t;ri 

-l 

1 

V( P ,e x=±1 )\Ql l (z)[z,-z}Q^-z)\o) = (pe)pJ-f 5 -f J ' d£e*M M ) 



for the XcO _meson: 



( X co(p)\Q l a (z)lz,-z}Q 3 p(-z)\o) = {P)fJ-f 5 f ld£e*M<fo&ri 



for the x cl -meson: 



for the L-meson: 



! 4 3 



fh k 



{h c (p,e\=±x) \Q(z)a flv [z,-z]Q(-z)\0) = {Phh) f3 a -^- ~^ 



for the Xc2-nieson: 



(xc2(p,tx=o)\Q l Q (z)lz-z}Q> p (~z)\o) = (p)pj-f 5 -f I dteVMfofan), 



Xc2(p,e x = ±1 )\Qi(z){z,-z]Q j /3 (-z)\o) = M x {pp) p J-f 5 -f 



The factor [z, — z], that makes matrix elements (|A1[) gauge invariant, is defined as 

[z, — z] = Pexp[ig f dx^A^(x)]. 
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In the above equations p is the charmonium momentum, x and x are the momentum fractions of quark and antiquark, 
£ = x — x, e M is the polarization tensor for the J/ip, Xd or he mesons and the vector p^ in relation (|A1|) is defined 
according to 

Pm = w 

where e M „ is the polarization tensor of the Xc2 meson. In practical applications it is useful to write the polarization 
of the \2 meson in terms of the polarization of two vector mesons. Thus, for instance, the polarization tensor of 
the transversely polarized \2 meson can be written as £x=±i = ( e A=±i ' e A=o + e x=o ' e \=±\)/^ ( e i/ e/ "' = •"■)' ^ 
we further contract the polarization tensor with lightlike four- vector z, to the leading twist accuracy we will get 
e Ml ' z v = e^_ ±1 (pz) I (y/2M X2 ) or p^ = e^_ ±1 /(v2-Mx 2 ). This form of the vector p can be used in the calculation with 
the leading twist accuracy. 

It is not difficult to show that the functions 4> v (0i </>^' T (£), an d 4>h(0 are £-even. The normalization 

condition for these functions is 

1 

J HOd( = 1. (A3) 
-l 

The functions <^£o(£), ^i(0) 0ft (£) anc ^ (0 are £-odd and normalized according to 

l 

J &(0<% = i. 

-1 
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